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ABSTRACT 

We analyze the X-ray properties of a near-infrared selected (K s < 20) sample of spectroscopically identified Ex¬ 
tremely Red Objects ( R-K s > 5) in a region of the Chandra Deep Field-South using the public 1 Ms observation. 
The 21 objects were classified, on the basis of deep VLT spectra, in two categories: 13 dusty star-forming galax¬ 
ies showing [Oil] emission, and 8 early-type galaxies with absorption features in their optical spectra. Only one 
emission line object has been individually detected; its very hard X-ray spectrum and the high intrinsic X-ray lu¬ 
minosity unambiguously reveal the presence of an obscured AGN. Stacking analysis of the remainder 12 emission 
line objects shows a significant detection with an average luminosity Lx ~ 8 x 10 4<) erg s -1 in the rest-frame 2-10 
keV band. The stacked counts of the 8 passive galaxies do not provide a positive detection. We briefly discuss the 
implications of the present results for the estimate of the Star Formation Rate (SFR) in emission line EROs. 

Subject headings: galaxies: active - galaxies: starburst - surveys - X-rays: galaxies 


1. INTRODUCTION 

Having the colors expected for high-redshift passive ellip¬ 
ticals, Extremely Red Objects (hereinafter EROs; R — K> 5) 
can be used as tracers of distant and old spheroids as a test for 
different cosmological models. The optical and near-infrared 
properties of EROs are also consistent with those expected for 
high-redshift dusty starbursts (e.g., HR 10, Cimatti et al. 1998; 
Smail et al. 1999) and Active Galactic Nuclei (AGN) reddened 
by dust (Pierre et al. 2001), for which the spectral energy dis¬ 
tribution at these wavelengths is dominated by the host galaxy 
starlight (e.g., Mainieri et al. 2002). 

The relative fraction of old ellipticals, dusty star-forming sys¬ 
tems and obscured AGN among optically selected ERO sam¬ 
ples is a key parameter in the study of the galaxy evolution and 
can constrain models which link the formation of massive el¬ 
liptical galaxies and the onset of AGN activity (Almaini et al. 
2002; Granato et al. 2001). A step forward in this direction 
has been recently achieved by the K20 survey (Cimatti et 
al. 2002a; Daddi et al. 2002). Deep VLT spectroscopy of a 
bright subsample at K s < 19.2 has shown that EROs are nearly 
equally populated by old, passively evolving systems and dusty 
star-forming galaxies over a similar range of redshift ( z m ean — 1 
for both the classes). 

Sensitive hard X-ray observations provide a unique tool to un¬ 
cover hidden AGN activity among the EROs population. In¬ 
deed, a sizeable fraction (10-30%) of hard X-ray sources re¬ 
cently discovered in deep Chandra and XMM surveys are as¬ 
sociated with EROs and the fraction depends on the limiting 
fluxes reached both in the optical and X-ray bands (Alexander 
et al. 2002; Mainieri et al. 2002; Crawford et al. 2001). The 
hard X-ray spectra indicate that X-ray bright EROs are associ¬ 
ated with high redshift, highly obscured AGNs. 


Moreover, thanks to the very deep limiting fluxes reached by 
deep Chandra surveys (CDFS, Giacconi et al. 2002, hereinafter 
G02; CDFN, Brandt et al. 2001a) it is also possible to in¬ 
vestigate via stacking analysis the average X-ray properties of 
fainter non-AGN sources (Brandt et al. 2001b; Hornschemeier 
et al. 2002; Nandra et al. 2002). 

On the other hand, the AGN fraction in optically selected EROs 
samples is still very uncertain and possibly lower than 15% 
(Brusa et al. 2002; Alexander et al. 2002), suggesting that the 
bulk of the EROs population is not related to active phenomena. 
In fact, the majority of optically selected EROs remains unde¬ 
tected even at the deepest limiting fluxes of the 1 Ms Chandra 
Deep Field North observation (Alexander et al. 2002). Stack¬ 
ing analysis of those objects not individually detected suggests 
a softer X-ray spectrum than that of AGN-powered EROs and 
an X-ray to optical flux ratio consistent with that expected from 
normal elliptical galaxies, star-forming systems or low lumi¬ 
nosity (L 2 ~\okeV < 10 42 erg s -1 ) nuclear activity at z — 1. 

Here we make use of the available spectroscopic information 
of a well-defined, magnitude limited ( K s < 20) subsample of 
EROs in the Chandra Deep Field-South area obtained as a part 
of the K2 0 survey (Cimatti et al. 2002a); the main aim of 
the present study is to investigate whether the X-ray emission 
properties of EROs are depending on the optical/infrared spec¬ 
tral properties. 

Throughout the paper, a cosmology with Hq = 70 km s -1 Mpc 1 , 
O m =0.3 and Q\=QJ is adopted. 

2. THE SAMPLE 

The ESO VLT K2 0 sample is one of the largest and most com¬ 
plete spectroscopic sample of galaxies with K s < 20 available to 
date; the main scientific goal of this project is the comparison of 
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the observed redshift distribution of a complete sample of about 
500 galaxies (87% complete to K s < 20, 98% if photometric 
redshifts are included) with the predictions of galaxy forma¬ 
tion models in order to obtain stringent clues on the formation 
and evolution of the present-day massive galaxies (Cimatti et 
al. 2002b). The K2 0 sample cover a sub-area (32.2 arcmin 2 ) 
of the Chandra Deep Field South (CDFS) and a 19.8 arcmin 2 
field around the quasar 0055-2659. 

Extremely Red Objects in the K2 0 sample were selected with 
the criterion of R-K s >5, and then classified on the basis of 
their optical spectra into two categories (Cimatti et al. 2002a): 
objects with an absorption line spectrum, consistent with that of 
early-type, passively evolving galaxies (hereinafter “old”) and 
objects showing [OIIJA3727 emission line with an observed 
Equivalent Width EW> 20 Aover a red continuum, consistent 
with that of dusty star-forming systems (hereinafter “dusty”). 
The K20 EROs sample in the CDFS area includes 48 objects at 
K s < 20. For the purposes of the present paper we consider only 
the 21 spectroscopically identified EROs which have been clas¬ 
sified in the two categories described above: 13 “dusty” and 8 
“old", distributed on a similar range of redshifts, z = 0.8 -f-1.6. 

3. X-RAY DATA ANALYSIS 

The Chandra Deep Field South (CDFS) survey is a combina¬ 
tion of 11 individual Chandra ACIS-I pointings observed in 
two years, from March 1999 to December 2000 (G02). The 
maximum effective exposure time is 942 ks, but it varies ac- 
cross the detector due to the CCD gaps and the rotation of the 
field-of-view in the different observations. In order to esti¬ 
mate the effective exposure time for each object, we have down¬ 
loaded from the archive, along with the combined photon event 
list, all the events corrisponding to the single pointings. We 
have calculated the exposure and instrument maps using stan¬ 
dard Cl AO 2.2.1 tools for each observation separately, and 
combined them to evaluate point-by-point the effective expo¬ 
sures. As the mirror vignetting is energy dependent, we cal¬ 
culated the exposure maps at two energies (1 keV and 4 keV), 
representative of the mean energy of the photons in the soft 
(0.5-2 keV) and hard (2-8 keV) bands. 

At first, we have searched for nearby (within 2") X-ray sources 
found by WAVDETECT in the K20 region, with a false-positive 
probability threshold of 10 -5 (Freeman et al. 2001). Only one 
source, a dusty ERO at z = 1.327, was individually detected; the 
optical position is almost coincident (A < 0.3”, within the ac¬ 
curacy of the K-band position) with the Chandra source CXO 
CDFS J033213.9-274526 in the G02 catalog. The X-ray and 
optical properties of this object are discussed in Sect. 4.1. 

In order to constrain average X-ray properties of the remain¬ 
ing individually undetected EROs, we have applied the stack¬ 
ing technique (Brandt et al. 2001b; Nandra et al. 2002) sepa¬ 
rately for the two classes of objects. The samples consist of 12 
“dusty” and 8 “old” EROs, with average redshifts of z = 1.053 
and z= 1 • 145, respectively. 

For the photometry we used a circular aperture with a radius 
of 2" centered at the positions of the selected EROs. We have 
also checked different aperture radii (1 -4”) and found that the 


adopted one is an optimal extraction radius for the counts esti¬ 
mate, in agreement with the Nandra et al. (2002) findings. The 
counts were stacked in the standard soft, hard and full bands 
(0.5-2 keV, 2-8 keV, and 0.5-8 keV), and in the 1-5 keV band 
which roughly corresponds to the rest-frame 2-10 keV band at 
the mean redshift of the objects. 

To verify whether the stacked net counts constitute a signifi¬ 
cant signal the local background has been estimated: in annular 
regions around each ERO centroid, as the average of eigth cir¬ 
cular regions (2 arcsec radius) for each ERO position and from 
a background map produced by the WAVDETECT algorithm. 
The different choiches do not significantly affect the results in 
none of the adopted energy bands, except the 2-8 keV range 
where the results are strongly affected by the increased level 
of the instrumental background. Extensive Monte Carlo simu¬ 
lations (up to 100,000 trials) have been also carried out shuf¬ 
fling 12(8) random positions for “dusty” and “old” EROs re¬ 
spectively, using the same photometry aperture (2 arcsec). The 
random positions were chosen to lie in “local background re¬ 
gions” to reproduce the actual background as close as possible. 
The resulting distributions are well approximated by Gaussians. 
The corresponding detection confidence level is in all the cases 
in good agreement with that computed from Poisson statistic. 
The results of the stacking analysis reported in Table 1 are re¬ 
ferred to the background estimated as the average of the eight 
circular regions; in Fig. 1 we show the summed images in the 
full (0.5-8 keV) band. An excess of counts above the expected 
background level is clearly present for the 12 “dusty” EROS. 
The signal is stronger in the full band (99.997% confidence 
level assuming a Poissonian distribution, corresponding to 
about 4.2a) and it is still present, although at a slightly lower 
confidence level, in the soft and 1-5 keV band, while it is not 
statistically significant in the hard band. Assuming a T = 2 
power law spectrum plus Galactic absorption (N w = 8 x 10 19 
cm -2 ), the 3a upper limit on the average band ratio (H/S<0.85, 
where H and S indicates the hard band and soft band counts, re¬ 
spectively), corresponds to an absorption column density lower 
than 10 22 cm -2 (see Fig. 4 in Alexander et al. 2002). 

Given that the X-ray spectrum of nearby starburst galaxies 
is well-represented by power-law models with slopes in the 
range 1.8-E2.5 and low intrinsic absorption (e.g., Ptak et al. 
1999; Ranalli et al. 2002a), we adopted I =2.1 to compute X- 
ray fluxes and luminosities. With this assumption, the stacked 
count rate in the 1-5 keV band corresponds to an average 2-10 
keV rest-frame luminosity of ~ 8 x 10 4 ° erg s -1 , at the mean 
redshift z=1.053. 

The X-ray emission from “old" EROs remains undetected 
in all the considered energy bands (see Table 1 for de¬ 
tails). Assuming a thermal emission model with kT=l 
keV and solar metallicity, fully consistent with the aver¬ 
age X-ray spectrum of nearby elliptical galaxies (e.g., Pel¬ 
legrini 1999), the 3a upper limit on the 0.5-2 keV lumi¬ 
nosity is Lx < 10 41 erg s -1 (see also Vignali et al. 2002). 



Stacked EROs in the CDF-S 
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Fig. 1.— Stacked full-band images of the “dusty” (left panel ) and 
“old” {right panel ) EROs. The images are 9x9 arcsec and have been 
smoothed with a gaussian of <r=1.5 pixel (approx 0.7 ,/ ) . The circles 
are centered on the stacking position and have a radius of 2". The de¬ 
tection significance of the summed counts is 4.2cr for the “dusty” and 
< 2cr for the “old”. The bright spot in the right panel at about A" 
from the stacked position, corresponds to an individually detected object. 


TABLE 1 

Stacking Results 



0.5-2 67 37.9 99.995% 10.4 1.3±0.45 7.4 ±2.6 

1- 5 83 61.4 99.6% 10.2 1.4±0.5 8.0 ±2.8 

2- 8 105 86.3 97.2% 10.1 <4.8 <27 



': Fluxes (in units of 10 17 erg cm 2 s and luminosity (in units of 10 40 erg 
s -1 ) are computed assuming a power-law spectrum with T = 2.1 and Galactic 
absorption for the “dusty” sample, and a thermal emission model with kT=l 
keV and solar metallicity for the “old” sample. Errors are at the 90% confi¬ 
dence level. 

*: For the “old” sample we report only the results in the 0.5-8 keV and 0.5-2 
keV bands. In the harder bands the results are always consistent, within 1 a, 
with the background fluctuations. 

4. RESULTS 

4.1. CXO CDFS J033213.9-274526: a type 2 AGN 

The only ERO individually detected in the 1 Ms exposure 
is the highest redshift, highest luminous (R=23.81, K s = 18.44) 
object in the “dusty” sample. It is a very hard X-ray source 
(F 2 - 8 fev=7.34 x 1(T 16 erg cm" 2 s“\ r = 1.4 and N H =Gal; G02) 
undetected in the soft band (Fo. 5 - 2 fev < 5.48 x 10 -17 erg cm -2 ) 
The observed band ratio (H/S> 7.3) implies an intrinsic X-ray 
column density larger than 4 x 10 23 cm -2 , for a classic AGN 
power-law spectrum with T = 1.8). The observed 2-8 keV lu¬ 
minosity is ~ 7 x 10 42 erg s -1 which corresponds to an unab¬ 
sorbed, rest-frame hard X-ray luminosity >4 x 10 43 erg s _1 , 
typical of Seyfert galaxies. At the high-energy, therefore, CXO 
CDFS J033213.9-274526 is naturally classified as a Type 2, ac¬ 
tive object. In the optical domain it shows only a faint, but 
significant, [OIIJA3727 emission feature. Unfortunately, at the 
expected [NeV] A3426 position a strong cosmic ray has affected 
the data. The lack of (broad) MgIIA2798 emission indicates 
that obscured activity is taking place also in the optical nucleus. 

1 SFR (M 0 yr 1 ) =2.0 XKT 40 U-lOkeV (erg s* 1 ) 


Given the small area sampled and the spectroscopic complete¬ 
ness, the present ERO sample is not suitable for a reliable esti¬ 
mate of the AGN fraction among the EROs population. How¬ 
ever, the present low detection rate of AGN among EROs (1/21) 
confirms that the bulk of the EROs population is not related to 
active phenomena. 

4.2. Dusty EROs 

The X-ray and optical data both suggest that starburst ac¬ 
tivity, rather than AGN, is powering “dusty” EROs. Indeed, 
following Cimatti et al. (2002a), we recomputed the average 
spectrum of the 12 EROs we used in the X-ray stacking anal¬ 
ysis. The average optical spectrum is in excellent agreement 
with that published in Cimatti et al. 2002a (see their Fig. 2), de¬ 
rived from the complete K20 EROs sample, except for the pres¬ 
ence of the AGN indicator [NeV]A3426; in the present “dusty” 
sample this feature can be safely ruled out (EW<1.5 A at the 
3er level), suggesting that the optical emission is dominated by 
star-forming systems. 

The inferred, average X-ray luminosity is around one order 
of magnitude larger than that of normal spiral galaxies (Mat- 
sumoto et al. 1997) and it is similar to that of the starburst 
galaxy M82 (e.g., Griffiths et al. 2000; Kaaret et al. 2001); 
the average X-ray spectrum (see Sect. 3) is consistent with 
that measured for nearby star-forming galaxies (e.g. Dahlem, 
Weaver & Heckman 1998; Ptak et al. 1999). 

We have also investigated whether the stacked signal could be 
due to a few individual sources just below the detection thresh¬ 
old. This is particulary relevant for our sample which may con¬ 
tain a few sub-threshold AGN responsible for the entire stacked 
signal. Five out of 12 “dusty” EROs are individually “detected” 
at the 95-99% confidence level. The X-ray luminosity for these 
“marginally detected” EROs never exceeds 6 x 10 41 erg s -1 in 
the full-band excluding the presence of a bright AGN among 
the stacked dusty EROs. Although the presence of a low- 
luminosity AGN (LLAGN) can not be completely ruled out on 
the basis of the present analysis, in the following we assume 
that the “dusty” EROs are mainly powered by starburst related 
processes. 

4.2.1. Star Formation Rate: the hard X-rays view 

There is an increasing evidence that the hard X-ray luminos¬ 
ity of star-forming objects can be used as a star formation rate 
(SFR) indicator (Ranalli, Comastri, & Setti 2002a,b; Nandra et 
al. 2002; Bauer et al. 2002). 

Using the Ranalli et al. relation 1 we have derived an average 
SFR for the EROs in the “dusty” sample of 5-44 M 0 yr 1 , 
taking into account the observed dispersion between the far- 
infrared and X-ray luminosities quoted by Ranalli et al., and 
the uncertainties on the hard X-ray luminosity. The average 
SFR derived from the the observed [Oil] luminosity (Kenni- 
cutt 1998) is Sf ; R(jOII])„/, s =2.7 M.. yr* 1 . The discrepancy 
between the optical and X-ray SFR estimates confirms that 
the optical emission lines are severely obscured by dust and 
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the optical value can be addressed as a lower limit. In order 
to provide a dereddened SFR estimate, the average spectrum 
was compared with template spectra of star forming galaxies 
and dusty systems (as extensively discussed in Cimatti et al. 
2002a); the best-fit estimates for the reddening are in the range 
E(B-V) = 0.5 4-0.8, which implies an intrinsic average SFR of 
the order 37-110 M 0 yr 1 according to the Kennicutt (1998) 
relation. Although the two estimates suffer of large uncertain¬ 
ties, mainly related to the scatter of the various SFR indicators, 
they can be reconciled if the lowest extinction correction (e.g., 
E(B-V) ~ 0.5, or even lower) is applied to correct the [Oil] 
luminosity. 

From the derived X-ray SFR it is possible to estimate the av¬ 
erage far-infrared (60-100 /im) and radio (1.4 GHz) emis¬ 
sion of the objects in our sample: Lfir = 10 10,4 4-10 11 ' 4 L 0 ; 
Li. 4 CHr = 2 4- 17.6 x 10 29 erg s -1 Hz -1 (Kennicutt 1998). The 
relatively low value of L/r /s would naturally explain the low 
detection rate of EROs in submillimeter continuum follow-up 
observations, typically sensitive to detect Ultra-Luminous In¬ 
frared Galaxies (ULIGs, L> 10 12 L 0 ) at z>l (e.g., Mohan et 
al. 2002). The “dusty” ERO population would then be rep¬ 
resentative of less luminous dusty star-forming systems, with 
respect to submillimeter-selected high-z galaxies, as already 
pointed out in Cimatti et al (2002a). 

The inferred radio luminosity corresponds to K-corrected 
(S„ oc v~ a , a = 0.7) radio fluxes in the range 4-38/iJy. The lack 
of deep radio observations in the Chandra Deep Field South 
area prevents to test these predictions; however, it is interesting 
to note that very deep 1.4 GHz VLA observations of a sam¬ 
ple of optically selected EROs (Smail et al. 2002) agree with 
our estimates (about 85% of the EROs in their sample is not 
detected in the radio band at fluxes higher than ~ 40 /tJy). 


5. SUMMARY 

Thanks to the deep limiting fluxes reached in the 1 Ms ex¬ 
posure of the Chandra Deep Field South and to the subarcsec 
positional accuracy of the X-ray satellite, it was possible to 
constrain via stacking analysis the high-energy properties of a 
well-defined, spectroscopically identified sample of EROs, be¬ 
longing to the K20 survey. 

The dichotomy in the optical classifications of EROs seems to 
be present also in the X-ray band: “dusty” objects are relatively 
bright X-ray sources, while “old” EROs are below the detection 
threshold. 

The only object (CXO CDFS J033213.9-274526) individually 
detected by Chandra is a luminous, highly absorbed AGN 
(N h > 4 x 10 23 cm -2 ) at z = 1.327. 

Using the hard X-ray luminosity as a SFR indicator, which is 
essentially free from intrinsic absorption, we have estimated a 
relatively low SFR (5-44 M 0 yr -1 ) with respect to the estimate 
based on the reddening-dependent [Oil] emission; our estimate 
is consistent with the results of recent far-infrared and radio ob¬ 
servations. 

Deep X-ray observations of large and complete samples of 
spectroscopically identified EROs could provide a powerful 
tool to estimate, via stacking analysis, the average star forma¬ 
tion rate in “dusty” systems and to better assess the AGN frac¬ 
tion among the ERO population. 
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